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Vertical Melting of a Stack of Membranes 
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Abstract. A stack of tensionless membranes with nonlinear curvature energy and vertical harmonic inter- 
action is studied. At low temperatures, the system forms a lamellar phase. At a critical temperature, the 
stack disorders vertically in a melting-like transition. 

PACS. 82. 65. Dp Thermodynamics of surfaces and interfaces - 68.35.Rh Phase transitions and critical 
phenomena 
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I ' 1 Introduction of each membrane to be less than the distance to its near- 
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• • i Under suitable conditions, lipid membranes in aqueous so- 

^ ' 

:<"^ ■ . , , , r i ii i • In this note, we investigate the question whether the 

K> , lution are known to form lamellar structures, cnaracteri- 

. , i i, , , r ,i i i, lamellar structure exists at all temperatures, or whether 

zed by a parallel arrangement ol the membranes alter- 

„ jTi, ™, . . r , thermal fluctuations can destroy the vertical order at some 

natmg with tnm layers ol water . 1 he existence ol sucn 

. . , ,1,1 r i critical temperature, 

structures is m contrast to the behavior ol a single ten- 
sionless membrane subject to thermal fluctuations, which 
is always in a disordered, crumpled phase, filling the em- 
bedding space completely 0. In a stack, this phase is sup- 
pressed by the steric repulsion between the membranes 
which prevents them from passing through each other jjj , 
thus constraining the amplitude of the height fluctuations 
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Such a transition does not take place in the simplest 
model of a stack proposed by Helfrich |jj , where the mem- 
branes possess only a linearized curvature energy and a 
harmonic repulsive term accounts for the vertical forces 
in the stack, approximating in a rough way the steric re- 
pulsion. In this purely Gaussian approximation, the the- 
ory is equivalent to de Gennes' theory of smectic-A liquid 
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In this paper we extend the simplest model by taking brane, and g m ,ij — % + diU m djU m the induced metric, 

into account higher-order terms of the curvature energy, withz,j = l,2,9i = d/dx,d2 = d/dyandg m = det[g m y]. 

and show that thermal fluctuations cause a finite stack The parameter tq is the surface tension of a single mem- 

of membranes to disorder vertically, a process which we brane, kq its bending rigidity, and Bq the compressibility 

may call vertical melting. A similar phase was previously of the stack. In Eq. ([!]), as in the following, the subscript 

described in a different context by Huse and Leibler Q, denotes bare quantities, whereas renormalized parameters 

and by Kleinert who related the molten phase to the will carry no subscript. 

directionally ordered phase of a nematic liquid crystal. We In the original Hclfrich model, the surface tension ro 

shall find that thermal fluctuations renormalize both the was not included because the membranes in the stack are 

coefficient of the curvature term and the vertical compress- tensionless. We have included ro in the energy (|l|) as an 

ibility of the stack. Analyzing the renormalization group infrared regulator which also serves to absorb ultravio- 

equations, we find, besides the Gaussian fixed point gov- let infinities. After carrying out the various integrals, the 

crning the low-temperature phase, a nontrivial fixed point renormalized physical tension will be set equal to zero, 

determining the critical exponents of the vertical melting For slow spatial variations, the discrete variable ml 

transition. may be replaced with a continuous one, and u(xj_, ml) —* 

it(x), where x = (xj_, z). In this limit, the energy (fil) re- 



duces to 

E = J " dz J d 2 x ± V5 [o-o + \K H 2 + \B Q {d z uf] . 

(2) 

Here we have introduced bulk versions of the surface ten- 
sion (To = ro/l and bending rigidity Kq = kq/1, and 
defined L\\ = Nl. In Ref. || the vertical gradient en- 
ergy (d z u) 2 was replaced by the normal gradient energy 
(N • Vu) 2 , which is physically more correct and has the 
advantage of being reparametrization invariant (see also 
Ref. §). In the following, we shall derive all results for 
both terms and analyze the difference between the two. 
(1) We study the stack perturbatively, starting with the 
Here, H m = diN m , is the mean curvature, where N m oc lamellar configuration, and expand the theory in the in- 



2 The Model 

We consider a generalization of the Helfrich model due to 
Janke and Kleinert ||, in which a multilayer system is 
made up of {N + 1) fluid membranes, parallel to the xy 
plane of a Cartesian coordinate system, separated a dis- 
tance I. If the vertical displacement of the rath membrane 
with respect to this reference plane is described by a func- 
tion u m (x) = w(xj_, ml), where x^ = (x, y), the energy of 
the stack reads: 



E = ^ J d 2 x^^/g^, 



r + ^Ko-ffm + -^T ( U m - «m-l) 2 



-diu m , —d2U m , 1) is the unit normal to the mth mem- verse parameter ao = 1/Kq, which is assumed to be small. 
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To keep our notation in conformity with the literature, we in the z-direction, the Fourier transform includes a sum 
use Kg throughout the text and resort to its inverse only (1/L\\) Yln=-N/2 over * ne discrete wavevector components 

when necessary. 27r 

u„ = —n. (4) 

L \\ 

We take into account the interlayer spacing I in a rough 
way by restricting the values of the discrete variables to 
\w n \ < tt/1, so that the summation index n lies in the 
interval — iiV < n < ^N. To one-loop order, the bare 
parameters are renormalized to 



3 Renormalization 

The phase transition to be described in this note is caused 
by a competition between the softening of the bending 
rigidity due to thermal fluctuations and the stack-preser- 
ving vertical elastic forces. To understand this competition 
we study how thermal fluctuations renormalize the param- 
eters of the theory. For this we expand the energy (0) up 



(To^<7o(l+/i), K -^K (l-%I 2 ), B -> B (l±il 2 ), 

(5) 



where 

k B T N/2 f d 2 q^ \—< jj2 - q±- f — q± 

to fourth order in the displacement field it(x), where it I\ = — — > / 7- .„ ° a „ ^ °° A (6) 



reads 



in , 

2, 



E= / dz/d^_L 







i a Q (d t u) 2 + ^K (dfu) 2 + \B {d z uf 



-—N/2' 

2 L W n= N/2 J W B0UJ " + ^ + K ^ ' 



We regularize the integrals in the ultraviolet by introdu- 
cing a sharp transverse wavevector cutoff A inversely pro- 
- K Q (d t u){d u){did.ju){dlu) portional to the lateral size a of a lipid molecule. Actually, 

± jBo(d z u) 2 (diu) 2 , (3) the divergent contributions to the above integrals are in- 

dependent of Bq. This is due to the discrete nature of the 

where the lower sign in the last term refers to the more 

stack. By restricting the values of the discrete wavevec- 

physical normal gradient energy (N- Vu) 2 . For zero surface 

tors uj n to account for the interlayer spacing, as explained 

tension, the lowest-order contribution to the energy due 

above, all divergences proportional to Bq are suppressed. 

to longitudinal displacements is of the usual elastic form 

The renormalization flow is obtained by integrating 

^Bo(d z u) 2 , while that due to transverse displacements is 

out transverse wavevectors in a momentum shell A/s < 

of higher order, viz. ^Kq{8 2 u) 2 . 

q± < A, and subsequently rescaling the coordinates. We 

The one-loop contributions can be calculated either by 

thus obtain, 

using Feynman diagrams, or by a derivative expansion, 

as in Ref. ||. Following standard procedure, we integrate ^ 1 ^ 2 Ly K ^ ^ U S 

out fluctuations with transverse wavevectors in a momen- Since these results are independent of 00, we can safely 

turn shell. Since the stack is periodic and of finite extent set the surface tension to zero, thus describing a stack of 
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tensionless membranes, characterized by the two remain- 
ing parameters, B and K. Under a rescaling xi — > xj_/s 
of the coordinates in the plane and z — > z/s Zc along the 
stack axis, the expansion parameter of perturbation the- 
ory a = l/K scales like a — > s~ Zc a, and the compressibil- 
ity scales like B — ► s 4 ~ Zc B. Here, z c allows for the pos- 
sibility of anisotropic scaling. From Eq. (g) with a = 0, 
it follows that z c = 2 in the Hclfrich model. Using the 
above results, one readily generates differential recursion 
relations to lowest nontrivial order 



da 
d In s 
dB 



3 k B T 



4tt U 



a 2 (N + l) 



(9) 



= (A-z c )B± 



1 k B T 



Ba(N + l). (10) 



dins 4-7T Lii 

Besides the Gaussian fixed point (a = 0, B = 0) which 
is stable in the infrared, the flow equations also admit a 
nontrivial, unstable fixed point at 

47T ill Z c , 



B = 0. 



(11) 



3 k B TN + l 

The latter implies the presence of a phase transition at a 
critical temperature 

4tt N lz c 



kv,T r 



(12) 



3 N + 1 a* 

The flow diagram corresponding to the above system 
of differential equations is shown in Fig. [I]. For T < T c , 
the inverse bending rigidity a flows to the Gaussian fixed 
point at the origin. In this low-temperature phase, the 
bending rigidity of the membranes increases with their 
lateral size, and thermal fluctuations are suppressed. This 
weak-coupling phase is the lamellar phase, where the ro- 
tational symmetry is spontaneously broken. For T > T c , 
on the other hand, a flows with increasing length scales 




Fig. 1. Flow diagram in the (a, _B)-plane. The diagram is plot- 
ted using the lower sign in the last term in Eq. (|ici|). 

away from the nontrivial fixed point at a* in the opposite 
direction. As a increases, the bending rigidity decreases 
and the membrane fluctuations become stronger. At the 
critical point a* , the stack disorders vertically and the sys- 
tem enters a strong-coupling disordered phase. Note that 
the critical temperature ( |l2| ) depends only weakly on the 
number (iV + 1) of membranes. 

The flow equations @) and ( |Io|) can be integrated ex- 
actly, yielding: 

* 4/* c -l±l/3 



B 



,±1/3 



(13) 



where c is an integration constant. For z c — 2, the expo- 
nent is equal to 1± 1/3. Explicitly, as we approach T c from 
below, B goes to zero as 



B ~ \T-T r \ l±1 / 3 . 



(14) 



The free energy density of the model in the lamellar 
phase can be calculated in the harmonic approximation, 
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as in Ref. [hi. For a finite stack, it reads 



/ = 



1 fc R T ( B 



16tt L 2 \K 



1/2 



- N(N + 2) 



(15) 



Thus, as T approaches T c from below, the free energy 
density behaves for z c = 2 like 

/HT-Tel 1 /^ 1 / 6 , (16) 

and the specific heat of the stack diverges as 

C~ |T-T c |- 3 / 2±1 / 6 . (17) 

Figure || shows a plot of the free energy density and spe- 
cific heat of the stack for the lower sign in Eqs. ( |To| ) and 

m 



c 



Fig. 2. Free energy density and specific heat of a stack of 
membranes, with the lower sign in Eqs. © and (0). 

We thus see that by using either the vertical gradient 
energy (d z u) 2 or the more physical normal gradient en- 
ergy (N • Vu) 2 , the qualitative behavior of the stack of 
membranes is not altered, but the critical exponents of 
the melting transition differ from each other. This is due 
to the fact that the normal gradient energy is zero for the 



in-plane flow of molecules inside the membranes. The ver- 
tical gradient energy (d z u) 2 , on the other hand, includes 
the energies of the tangential flow. Incompressibility ef- 
fects have been shown by David [|l0| and by Kleinert [|ll] 
to be irrelevant for the renormalization of a single mem- 
brane. Our result implies that this is not the case for a 
stack of membranes. 

The properties of a single membrane are obtained by 
letting N — > and z c — > 0. In particular, the flow equa- 
tion @ of the bending rigidity reduces in this limit to 
the known result . It has no fixed point other than the 
trivial one, which is unstable now, implying that a sin- 
gle membrane is always in the crumpled strong-coupling 
phase. 

4 Weak-coupling ordered phase 

We next wish to study the low-temperature phase of the 
stack in more detail, where the coupling constant a is 
weak. As in smectic-A liquid crystals, the theory (^5j) can 
here safely be approximated by taking into account only 
the quadratic terms. An important characteristic is the 
vertical fluctuation width or roughness t of a membrane 
in the stack defined by the mean square height fluctuation 
as £ 2 = {u 2 ) |12] ]. In the harmonic approximation, it is 
given by the one-loop integral 

N/2 

JT> (I 8 ) 



2x _ k B T y-v 



L 



-—N/2 



d 2 9J 



1 



1/L± (2tt) 2 Blu 2 + Kqj 



where, in the absence of a surface tension, the largest 
wavelength is equal to the inverse lateral size 1/L± of 
the membranes. As first observed by Peierls and Landau, 
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the mean square fluctuations diverge in the infrared. They going from the limit of a single, tensionless 2-dimensional 
thus destroy the long-range positional order in the layered (D = 2) membrane to the opposite limit of an infinite 
system at any finite temperature. More specifically, one stack, we go from the lower critical dimension of the for- 



finds jT2 13 1 mer to the upper critical dimension of the latter 

n2 k B T 



2n 



(19) Another characteristic of the weak-coupling phase of 
low-temperature is the behavior of the structure factor 



KL\\ VBK 

The first contribution, only present in a finite stack, stems 
from the n = term in the sum in Eq. (|l8|). It corresponds 
to a super soft mode, where the membranes undulate 

coherently with constant interlayer distance. The second where q is parallel to the z-axis, with \q \ = 2n/l. This 



S n (x) = (exp{zn go [u(x) - u(0)]}), (21) 



contribution, on the other hand, is also present when the correlation function can be directly observed in X-ray scat- 
stack is infinite. This contribution increases slowly with tering experiments, where the fluctuation spectrum is ex- 
the lateral size. pressed as half-widths at half-maximum of the anomalous 
Before proceeding, let us pause for a moment and con- Bragg peaks. As in smectic-A liquid crystals the 
sider the roughness of the two limiting cases of our theory: Fourier transform of the structure factor has algebraic sin- 
a single, tensionless membrane and an infinite, continuum gularities at q z = nqo: 

stack of such membranes. In this way, we find 

d D q±k B T S n (0, qz )~(q z -nq Q )- 2+n2 \ S(q±, 0) ~ q^ 2 ^, 

TTTTn — T' (single membrane) 

(27r)^ nqj_ 

dq z d D q± k B T ...... . ( 22 ) 



(2.)^ Bql+Kqi' (hlfinite StaCk) ' 

(20) 



with exponent 77. In the harmonic approximation, rj can 
be calculated from Eq. (|2l]) and turns out to be the same 
as for an infinite stack fll4| 



where instead of a 2-dimensional membrane we consider 
a Z?-dimensional object. It follows that for D > D u = 

4 — z c , the roughness is finite in the infrared, indicating k B T q 2 

that D u = 4 — z c is the upper critical dimension. Recall 



8tt JBK 



(23) 

that for a single membrane z c = 0, while for an infinite The a i ge braic singularities in @ reflect the quasi-long- 
stack z c = 2. To determine the lower critical dimension, we range periodic order along the stack axis. As for smectic-A 
consider the mean square normal, or omentational fluctu- liquid crysta ls, the exponent r) is temperature indepen- 
ations ((<9 ± u) 2 ). This results in an additional factor of q\ dent. This can be seen by remembering that by simple 
in the numerator of the integrands in (pOj). The resulting scaling arguments 
expressions are finite in the infrared for D > D\ = 2 — z c , 

p . (fc B T) 2 I 

identifying D\ as the lower critical dimension. Hence, in K n _ w \4 ' ^ ' 
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for a stack of membranes of rigidity k and thickness w. 
Specifically, 

,~o-?r. (25, 

For smectic-A liquid crystals, this expression was con- 
firmed experimentally |lq ]. 

As the temperature increases, thermal fluctuations be- 
come stronger, and eventually overcome the vertical forces 
of the B-term leading to a vertical disordering of the stack. 
In the strong-coupling disordered phase, the Gaussian ap- 
proximation used in this section to investigate the weak- 
coupling ordered phase breaks down. To study this phase, 
a nonperturbative method is required. 

5 Conclusions 

We have shown that a stack of membranes with nonlin- 
ear curvature energy melts at some critical temperature 
T c . In the weak-coupling low-temperature phase the sys- 
tem forms a periodic array of well-defined surfaces. There 
is long-range orientational order in the planes, and quasi- 
long-range positional order along the stack axis. This phase 
can be accurately described by the harmonic approxima- 
tion of the Helfrich model, which coincides with the de 
Gennes' theory of smectic-A liquid crystals. 

Upon approaching T c , the stack melts and the lamellar 
phase goes over into a strong-coupling disordered phase. 
This phase cannot be described by the harmonic approx- 
imation. Its properties will be investigated separately in 
a nonperturbative framework, in the limit of infinite em- 
bedding space dimension. 
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